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A scaling theory of replica symmetry breaking (RSB) in the SK-model is presented in the frame- 
work of critical phenomena for the scaling regime of small inverse RSB-orders 1/k, small tempera- 
tures T, and small magnetic fields H . We employ the pseudo-dynamical picture (Phys. Rev. Lett. 
98, 127201 (2007)) with two critical points CV1 and CV2 where separated temperature- and mag- 
netic field-scaling is obtained near the order function's pseudo-dynamical limits Hm a ^ oa q(a) — 1 
and lim a _ r0 q(a) = at (T = 0, H = 0). An unconventional scaling hypothesis for the free energy is 
given, modeling this separated scaling in accordance with detailed numerical self-consistent solutions 
for up to 200 orders of RSB. Divergent correlation-lengths £cvi(T) ~ T~ v t and &V2(H) ~ H~ v h 
describe the RSB-criticality as a long-range correlation effect occurring on the pseudo-lattice of 
RSB-orders. Rational- valued exponents v T = 3/5 and u H = 2/3 are concluded with high precision 
from high-order RSB scaling (in analogy with finite size scaling) and using a new fixed point extrap- 
olation method. Power laws, scaling relations, and scaling functions are analyzed. Near CV1, the 
non- equilibrium susceptibility is found to decay like Xi = k -5 ^ 3 fi(T/ k~ 5 / 3 ), the T = 0-entropy like 
S ~ x?i while (T-normalized) Parisi box sizes diverge like dj = K 5//3 / ai (T//t~ 5//3 ), with /i(C) ~ C 
and / 0j (C) ~ 1/C f° r C °°j /(0) finite. Near CV2, where the magnetic field H controls the critical 
behavior (while temperature is irrelevant), a power law H 2 ^ 3 is retrieved for plateau-height of the 
order function q(a) according to q p i(H) = ft -1 f p i{H 2 ^ s /k _1 ) with / p i(C)lc->oo ~ C an d f P i(0) finite. 
The order function q(a) links CV1 with CV2 and is obtained as a fixed point function q*(a*) of RSB- 
flow, in agreement with integrated fixed-point energy and susceptibility distributions. Similarities 
with directed polymers in 1 + 1 dimensions, with d = 1 solution and Flory-Imry-Ma type solutions 
of the KPZ-equation are discussed. 

PACS numbers: 75.10.Nr,75.40.-s,89.75.Da 



I. INTRODUCTION 

The far-reaching usefulness of spin glass theoriesi^ 
and of its key structural elements such as frustration, 
disorder, hierarchical order, ultrametricity, complexity, 
freezing transition, is witnessed well by applications en- 
tering even life-sciences and trans-disciplinary research 
fields. Physical models, where these key structures ac- 
quired a specific mathematical meaning, find very broad 
applications beyond their origin in frustrated magnetism. 
Let us mention, apart from fields like neural networks, 
computer science, and econophysics, the fascinating so- 
ciological applications to opinion- and group dynamics^, 
biological applications to RNA-folding^i^ including the 
quantum chromodynamical analogy and random matrix 
theory^. It seems natural to search for universal features 
of unifying models both in the general sense and in the 
precise meaning of the renormalization groups. 

The 3SAT optimization problem and its close rela- 
tion with the T = Sherrington-Kirkpatrick models or 
RNA-folding in biophysics^, where glass transitions exist 
within the secondary RNA-structure^ i 10 ' 11 , provide ex- 
amples where even the zero temperature limit is cither 
exact or close to the realistic situation. In physics, spin 
glass phases are usually confined to a low temperature 
regime and some applications are rather remote from it. 
Yet, knowing the ground state structure remains impor- 
tant. For one of the most fertile standard models, the 



Sherrington-Kirkpatrick models (SK-model), the hier- 
archical ground state structure is meanwhile confirmed 13 
as predicted by Parisi a long time ago 14 i 15 . Explicit 
analytic solutions on the other hand or meaningful ap- 
proximations are still required. They may lead to im- 
proved understanding and could be potentially fruitful 
for progress in more complicated (non-mean-field finite- 
range, quantum-) models. 

The attempt to link the SK-model behavior deep in- 
side its ordered (spin glass) phase with the theory of crit- 
ical phenomena may appear unmotivated at first sight, 
since the infinite-ranged spin interaction suggests 'only' 
mean-field behavior. However the SK-model is not simple 
below its mean-field transition. Its replica theory^ al- 
lows to imagine how the Ising Hamiltonian with infinite- 
ranged random interaction can become potentially criti- 
cal, when it is dressed-up with the hierarchical order pa- 
rameter structure in the replica symmetry broken (RSB) 
phase 3 -. We shall argue in this paper that, as the number 
of tree levels of this hierarchical structure grow to infin- 
ity, a particular correlation length between the tree levels 
can be defined which diverges as T — > 0. It allows to de- 
scribe critical behavior due to the accumulated effect of 
ever finer structures at the highest tree level. This prop- 
erty specifies a kind of universality class, which helps to 
compare with similar behavior in different physical sys- 
tems and in other scientific areas like biology, sociology, 
(mathematical) psychology as well, where evidently frus- 
trated random (and in cases range-free) interactions are 



2 



important. 

Nonanalytic power laws (with rational exponents) for 
the SK-model had been discussed in many different re- 
spects, as for example the finite-size cutoff (or finite spin 
number) dependenc e 17 ' 18 ! 19 . One may also mention the 
exponent of the Almeida Thouless linei 6 -. However a link 
to specific critical points was not made. 

In the present paper, we shall report progress in under- 
standing replica symmetry breaking in the Sherrington- 
Kirkpatrick modeU 2 - as a critical phenomenon; this refers 
to scaling behavior on one hand and to (numerically de- 
termined) fixed point functions under RSB-order flow 
k — ► oo on the other. Nonanalytic scaling behavior is 
described as a function of the inverse RSB-order decreas- 
ing to zero either together with temperature T — > 0, or 
together with the external magnetic field. Temperature- 
and field-scaling are well separated and reside in opposite 
limits of a pseudo-dynamic variable 1/a (see Ref|20|), as 
sketched by Figfjp. 

Critical phenomena arc in general categorized by uni- 
versality classes and described by criteria like global sym- 
metries. Certain details (on shorter range) become ir- 
relevant and suppressed in the regime of divergent cor- 
relation lengths. In the early years of the develop- 
ment of phase transition theory and critical phenomena, 
Kadanoff's initializing ideas of universality and rescal- 
ing, Stanley's scaling theory, and Wilson's rcnormaliza- 
tion group led to the modern understanding of critical 
behavior—. In recent years the functional rcnormaliza- 
tion group was advocated to understand better disorder- 
related criticality22. 

Freezing transitions into spin glass phases were an- 
alyzed in renormalization framework too, the ordered 
phase itself remained however mysterious, in particular 
for the non-mean-field models. In a famous work on 
scaling in spin glasses D. Fisher and Sompolinsky 2 ^ ex- 
plained the complications of mean field models (or mean 
field regimes of finite-range spin glasses above d = 6 and 
d — 8) and the multiple violations of scaling relations. In 
particular they mentioned the violation of temperature- 
versus magnetic field scaling within the ordered phase. 
In a different manner, we re-encounter this problem and 
explain a certain decoupling of field- from temperature- 
scaling by the presence of two different critical points of 
RSB in the low temperature limit. 

Crucial questions like the relationship between Parisi's 
RSB and the Fishcr-Huse droplet theory 2 ^ of the ordered 
phase of real spin glasses (or their reconciliation) became 
- since a long time but perhaps currently with more good 
reason - a field of intensive researc h 25 ' 26 . Since droplet 
theory is interpreted to govern the ordered phase by a 
T = fixed point, it appears very desirable to understand 
RSB as a T = fixed point theory too. Attempts have 
been independently undertaken by several authors and 
also in different fields of application, as the examples in 
Re^HHHHm! show. 

The latter point is elaborated in the present article. 
Despite the mean field character of the SK-model, RSB 



introduces apparently nonanalytic critical behavior of 
one-dimensional type (the unbroken replica-symmetric 
solution does not show any of these phenomena) together 
with special diverging correlation lengths. The challenge 
to handle RSB-effects correctly and to make the SK- 
solution a fruitful basis for real physical applications led 
us to a scaling theory intimately based on extreme high 
order numerical results. 

In previous publication s 20 ! 27 we reported the existence 
of two critical points and of discrete spectra which sur- 
vived in the limit of infinite replica symmetry breaking 
(oo-RSB) for the SK-model at T = 0, perhaps surpris- 
ing since the oo-RSB limit is in generally known only 
as the 'continuum limit'. Indeed, a continuum scaling 
theory, dealing with the T — ► 0— limit at K = oo) was 
published by Pankov^ recently. Its role and limitation 
to the temperature-controlled critical point CV\ has been 
addressed in our previous publication^ - together with a 
comparison of our work with the much older so-called 
PaT-scaling^i. In the present article we do neither use 
Pankov- nor PaT-scaling, but construct a different scal- 
ing approach, which includes RSB-order-scaling, and is 
exclusively guided by the theory of critical phenomena. 
In accordance with previous (naive) functional renormal- 
ization group arguments^ 7 - we analyze the approach to 
full RSB formation (n — > oo) not only at T = but 
also in the (H, T)-plane for small values of temperature 
T and magnetic field H and, of course, as a function of 
RSB-order. (neither real space nor real-time space are 
involved as a consequence of the SK- model's nature). 

We suppose here that RSB orders, counted by inte- 
gers 1,2, ...,k, can be viewed as equidistant sites form- 
ing a pseudo-lattice. In analogy with a real-space lat- 
tice, which needs to be infinitely large in order to allow 
diverging correlation lengths and hence support critical 
phenomena, the pseudo-length cutoff k must be sent to 
infinity. The known fact that increasingly high orders of 
RSB are needed (for good approximations) as the tem- 
perature decreases towards zero implies the role of T as 
an effective cutoff of nonanalytic behavior in the RSB- 
limit (T playing the role of a symmetry breaking relevant 
perturbation in standard critical phenomena). Thus it 
also inherits the idea of scaling RSB-order k with tem- 
perature T. Conversely, a maximum RSB order k serves 
as a cutoff of criticality. A speciality of RSB is that 
it appears in the shape of a pseudo-dynamical critical 
phenomeno n 20 ! 27 , which recalls the celebrated dynamical 
representation of Sompolinsky^ 2 -. A technically impor- 
tant difference however being the absence of a stochastic 
field, which we reserve for more complicated couplings to 
faster degrees of freedom^. 

A scaling theory, near T = in particular, is im- 
portant for several different reasons. First, it expresses 
the numerically determined features of the SK-model 
in a universal form, which helps to identify model- 
independent features and places the SK-model and its 
RSB in a wider context. Let us mention that directed 
polymers (or for example the queuing transition and 
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the totally asymmetric exclusion proces o i 35 , or certain 
partial differential equations, involve rational exponents 
as multiples of 1/3 too). The scaling theory also puts 
constraints on the shape of an effective field theory, ft 
has the virtue of isolating critical features which must 
be represented correctly by an effective theory that 
simplifies the SK-modcl^ The simpler theory should 
allow to control generalizations to finite range or other 
complications. The scaling theory offers also a special 
look on eventual scenarios of an RSB breakdown, as 
it may occur due to finite range interactions. The 
collapse of finite T c below a lower critical dimension, 
will eventually combine RSB-criticality with the freezing 
transition as T c = 0. 

The paper is organized as follows: 

Section [TT] and IIIII describe the basic elements of the 
present scaling theory. The spaces spanned by the scaling 
variables at both critical points are described in section 
HI] In section IIIII a correlation length is introduced on 
the pseudo-lattice of RSB-ordcrs (to our best knowledge 
for the first time) and, anticipating the self-consistent 
numerical results of the following Sections IIVIIVIIII (for 
details see Refl30|). the role of finite temperatures (or 
finite magnetic fields) as soft cut-offs of the divergence of 
this correlation length is explained. 

Section IIVI demonstrates how the order parameter 
function q(a) can be regarded and obtained as a fixed 
point function q* (a*) under RSB-flow K — > oo. 

Section |V] includes and combines finite temperature 
scaling near critical point CV1 with RSB-order scaling. 
Scaling functions are obtained, which fit the detailed data 
of 200 RSB-orders, and explain the non-commuting sin- 
gular limit k^oCjT^O. In a similar way. I VII includes 
magnetic field scaling near CV2. 

In Section IVIII we present unconventional scaling- 
contributions to the free energy, to the entropy, and in- 
ternal energy, which are compatible with the numerical 
self-consistent solutions. 

In Section [Villi the ground state energy distribution is 
given as a function of pseudo-time and also shown as a 
function of the (normalized) Parisi levels l/n such that 
the flow towards a energy-per-level fix point function re- 
sults as the RSB-order tends to infinity. 

In Section HXl we finally consider pseudo-dynamic scal- 
ing of the order function q(a) in the vicinity of both crit- 
ical points before concluding with details of q(a) as re- 
vealed by its derivatives in Section [Xj 



II. THE SCALING SCENARIO 

We introduce the (RSB-)scaling idea by viewing the 
formation of full RSB as a critical phenomenon with two 
critical points in the pseudo-dynamic limits a = and 
a = oo at T = 0, H = 0. We do not a priori impose a 
relationship between the two critical points, but consider 



the pseudo-dynamical crossover between them by means 
of the order function q(a) on < a < oo. Fig[T]illustrates 
the relative position of the two critical points and the sets 
of scaling variables near these points. 

In particular one may notice that the dynamical vari- 
able 1/a and the RSB-order k define a (l + l)-dimensional 
analogy of problems with one time- and one real space 
dimension. Since the free energy or internal energy are 
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FIG. 1: Scaling regimes near the two critical points of the 
SK-model at zero temperature, zero field, and infinite replica 
symmetry breaking. The critical points are separated by the 
full range of continuously distributed pseudo-times 1/a. 

integrals over all pseudo-times, as for example given be- 
low in Eqs. (f2~If . (|26|) . we do not start from a single scaling 
hypothesis for the free energy F. Instead we construct 
the scaling hypothesis for each of the two different scal- 
ing contributions, originating in these separated critical 
points and. As FigQ] shows, a different set of scaling 
variables should be used in order to match the numerical 
results. 

It is remarkable that temperature- and magnetic field- 
scaling become decoupled, because they belong to differ- 
ent scaling regimes. Scaling with the respect to the order 
k of RSB measures the approach of the equilibrium so- 
lution at k — oo (full RSB) and therefore can be viewed 
as a kind of non-equilibrium dynamics (in the sense that 
each finite order is unstable towards higher RSB-orders). 
Thus an element of dynamic scaling is contained. Using 
the pseudo-time 1/a as an additional scaling variable, 
we analyze the order function q(a; T, H) and its pseudo- 
dynamic scaling behavior. A dynamic crossover between 
the two critical points CV1 and CP2 is then described by 
means of q(a). Moreover, the order function is evaluated 
as a fixed point function of the RSB-flow letting k — > oo. 

The present scaling theory is then fitted to high preci- 
sion numerical data, which were obtained recently for the 
Sherrington-Kirkpatrick model given by the Hamiltonian 

H = ^2 Jij s iSj - H 2J Si 
i<j i 

with quenched, infinite-ranged, and Gaussian-distributed 
random couplings Jy (with variance J 2 /N) between clas- 
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sical spins s, = ±1. The method was described in Refl30l 
and will not be described again in this article. It al- 
lowed not only to go beyond earlier high-order studies^, 
but contained also new analytical elements. As a conse- 
quence we are able to predict the values of critical expo- 
nents, evaluate amplitudes, calculate analytical models 
of various scaling functions including cases with very sin- 
gular crossover. 

The numerical material includes the self-consistent 
solutions in all orders of RSB up to 

i) the current maximum of ft = 200 RSB at T = and 
H = 0, 

ii) 50 orders for a dense grid of finite temperatures in 
the range < T < 0.3 for H = 0, and 

iii) 20 orders of RSB for a dense grid of finite magnetic 
fields < H < 0.5 at zero temperature. 

We note that all energies are given in units of J. 



III. DIVERGENT CORRELATION LENGTH ON 
THE ONE-DIMENSIONAL PSEUDO-LATTICE 
OF RSB-ORDERS 

The high order self-consistent solutions led us to con- 
sider a pseudo-lattice of RSB-orders with unit lattice con- 
stant. The maximal order k, for which self-consistent re- 
sults have been obtained numerically, can be viewed as a 
(sharp) pseudo-length cutoff. In analogy with the well- 
known finite size scaling of critical phenomena, one may 
consider scaling by varying this finite maximum RSB- 
order k. Naturally this defines a one-dimensional prob- 
lem without translational invariance though. Moreover 
it is known that finite temperatures or finite fields serve 
each as a soft-cutoff for the maximal order of RSB, which 
is needed to obtain good approximations: the higher the 
temperature the less RSB orders are needed to obtain a 
certain quality of results^. 

In other words, higher orders become uncorrelated be- 
cause they only have weak and/or negligible effects. In 
this sense the correlation length of different RSB or- 
ders becomes cutoff by finite T or H . Anticipating our 
results below, this definition of the correlation lengths 
£ K (a,T, H) shows power law divergences with rational- 
valued exponents given by 

to>x = t*(a = oo,H = 0,T)~T- v *, v T = 3/5 (1) 

&v2 = Ua = Q,H,T = 0)~H-'"*, v H = 2/3 (2) 

In the chapters below, we shall find scaling functions 
of the form /(«/£«)■ Apparently the correlation length 
exponent vh cannot be simply related to vt, since in the 
vicinity of one of the two critical points one finds either T 
and no //-dependence (CP1) or vice versa (CP2) (unlike 
conventional scaling where T ~ H 1 /^ 5 ). The power laws 
(1) and (2) should also hold for a ^> £ and a <C £o 
respectively, where £o was introduced in RefH?] as a finite 



characteristic length (« 1.13) which sets the scale for 
crossover from almost linear regime, q ~ a, through a 
maximal-curvature crossover (a ~ £o) to 1 — q m I /a 2 
behavior of the order function q(a, T = 0, H = 0). 

Let us note some similarities with critical behavior of 
other systems. Very remarkable appears the fact that 
the exponent vt = §, describing the RSB-correlation- 
lcngth divergence at CP1, coincides with the value given 
by Garel and Orland for the variational domain-wall so- 
lution of (1 + 1)— dimensional directed polymers^. 

For general d < 2 these authors reported v = 3/(d+4), 
while a second solution called domain-solution yields 
v = 1/(4 — d). Hence, in (1 + l)-dimcnsional case of 
directed polymers, these two solutions give v = 3/5 and 
v = 1/3 respectively. If H 2 in our case could be scaled 
like T then our second correlation length exponent 
vh would agree with the domain solution for directed 
polymers (DP). Moreover, the roughness-exponent of 
the DP in 1 + 1 dimension assumes precisely this value 
2/3. Of course, these are only hints, and universality 
classes can coincide accidentally at integer dimensions. 
It appears however very tempting to analyze whether 
our two critical points CV1 and CV2, distinguished by 
the opposite pseudo-dynamic limits a = oo and a = 
respectively, find their counterparts in those domain- 
and domain- wall solutions of the directed polymers. 
Whether such a formal analogy holds (despite the 
infinite-range spin glass interaction), should probably be 
decided by comparison or mapping of the corresponding 
field theories. Searching for eventual relations in higher 
dimensions for finite range spin glass interactions is also 
an exciting question. We shall come back to related 
questions and pseudo-dynamic scaling in section HXl 



IV. THE T = ORDER FUNCTION AS A FIXED 
POINT FUNCTION q*(a*) IN THE RSB LIMIT 

(ft = oo) 

The idea of finding the T = order function as a 
fixed point function in the RSB-limit arose from renor- 
malization group arguments as designed in Ref. l27l . It 
reemerges now in a literally obvious way when we plot 
in Figure [5] the whole set of numerical self-consistent so- 
lutions {a\ K , (and for I = 1...K of all 
evaluated RSB-orders k = 1,2, ...,200. These data be- 
come dense for large enough k and approach the desired 
order function q(a) in the limit n — > oo, which can be 
viewed as a fixed point function q*(a*). 

The unusual form q* (a*) can be justified as follows: the 
parameters and qf^ can be viewed as functions of 
the continuous variable 1/k — > £ in the n — > oo limit. The 
fixed point solutions a* (£) and q* (£) can be combined by 
eliminating the variable £, which results in the special 
form q*(a*) where the variable itself is made up from 
continuously distributed fixed points. In the following we 
use q(a) and q*{a*) synonymously and distinguish them 
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FIG. 2: Main Figure and Insert present an (overlapping) wide-scale overview of all (T = 0, H = 0) selfconsistent solutions 
{a-l^ , an d {a; , <Zj }, where / = l.../t for each RSB-order «: = 1,2, ...,200. One color cycle from k = 1 (red, remote 

points) to the highly dense regime « = 200 (red, almost continuous curve) demonstrates how the true order function q(a) 
is approached by the almost continuous g; + i(a;) from below and by qi(ai) from above. The Insert shows the corresponding 
{l/a,g}-data and the function q(l/a) they converge too. The shape near r = - = 0, (CV1) is almost Gaussian (dashed black 
curve: exp(— 0.41/a 2 )) but crosses over for moderate and small 1/a to strongly non-Gaussian and a tail behavior q(r) ~ 1/r. 



only if necessary. 

At finite large orders one may define interpolating 
functions qi+i(ai) — * q<(a) and qi{a{) — > g>(a) which 
yield lower and upper bounds for the exact solution 
q(a) = q*(a*) at each value of a. Figure [2] illustrates 
that this channel between lower and upper bound be- 
comes extremely small for high orders k = O(10 2 ). An 
illustration of the exact q(a) being confined within such 
a channel as the RSB order k increases towards infinity, 
is provided in a more detailed way by zooming different 
regions of crossover between CV1 and CV2 in Figures EHU 

Figure [5] moreover shows deviations of q(l/a) from 
Gaussian behavior, which is a good approximation for 
small 1/a. In both representations q(a) and q(l/a) it 
illustrates the existence of special lines which terminate 
obviously in fixed points - in fact there is a hierarchy of 
fixed points lying dense on the interval < a < oo. We 
shall make explicit use of these fixed points below. 

Indeed, 200 calculated orders of RSB for T = al- 
ready yield an almost continuous function qi ^ a '+^'+ 1 ^ 

which finally turns into q(a) in the RSB-limit. Our pre- 
viously published analytical model function^ satisfies al- 
most perfectly this constraint; as mentioned in Ref. 



it turned out that a small 'mass' function w(a) in 



a ( t 2 \ , s 

■sja 2 + w(a) V a 2 + w{a)J 



models even the full crossover regime. The function w(a) 
tends to a small constant w(0) 0.067. This cutoff 
(of in general nonanalytic behavior for arbitrary param- 
eters a, 7) guarantees a strictly linear q(a) relation in 
accordance with our high order data. In the crossover 
regime between these two dynamic critical points, w(a) 
can be modeled to depress the maximum error of q(a) 
below O(10 -4 ) at each pseudo-time. A unique choice 
of w(a) is not yet found, but excellent fits are obtained 
with w(a) monotonically decreasing from w(0) ps 0.067 
to w(oo) = 0. Using the high order data we have thus 
been able to improve the analytic approximation of the 
T = order function q(a). 
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FIG. 3: Self-consistent {a, q}-data in the small-a regime 
are shown for all k = 1,2, ..200 RSB orders, together with 
an analytical model function q(a) (dashed grey line). So- 
lutions {cii iK , qi +lyK } (above the line) and {a; ifc , qi, K } (be- 
low) approach fixed points under constrained RSB-flow obey- 
ing k = I level number. Fixed point examples for 
k — l(m + l)/m — k/m — 1, 23 > m > 6 and k integer 
are shown (red dots) including the fixed point (0, 0) (RSB- 
flow along k = I + k). The analytical model function for q(a) 
matches well all fixed points. 



A. Fixed points calculated from the RSB-flow 
towards k = do 

The full set of self-consistent solutions for order pa- 
rameters qi and (T-normalizcd) Parisi box sizes ai = 
mi(T) /T\t = can be described by matrix elements 

Pi,k = {ai, K = a\ K \qi >K = q^} labeled by RSB-ordcr 
k and level number I. Since the number of ^-parameters 
exceeds by one the number of a/-parameters (in each 
order of RSB), a second complementary set of matrix 
elements pi iK = {a/. K , qi+i, K } should also be taken into 
account. These points pi tK and pi tK arc displayed in 
the Figures 2 — 6 and observed to approach the exact 
q{a) = q*(a*) along characteristic lines given below by 
Eq.(j4]) as k — > oo (jp from above and p form below q(a) 
since qi+i, K < qi, K ). 

The set of all RSB-solutions up to a maximum order 
k , as plotted in Fig[5]with a cutoff at n = 200, is then de- 
scribed by two triangular matrices with entries {et/. K , qi_ K } 
(or with {a; iK , <7; +l K }); the level-numbers I run from 1 to 
k for each RSB-order k. 

Along infinitely many lines in (I, K)-space - the leading 
ones are very clearly visible in Figures 13141 (and shown as 
calculated in Figs. 5, 6) - we observe very smooth behavior 
of slowly changing parameters (a; iK , qi >K ) which allow low 
order Pade-approximants to match these data and to join 
in fixed points p* of the order function curve for k = 
oo. A special case is the origin where the fixed point is 
obtained with the extreme accuracy of O(10 -13 )). 

Typical examples of such characteristic lines in (I, k)- 
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FIG. 4: Two examples of alignment of solutions {a, q} 
under constrained RSB-flow with n(l) (or 1(h)) towards 
two fixed points (big dark- red dots): (top) v*m=~i ,n=6 = 
{a*,q*) m=7 , n =e = (0.3820,0.2976) and (bottom) p* m=a , n=2 = 
(a*,q*)3,2 = (0.9133,0.6559). Pade-line intersections, shown 
for 6 in Fig[5]and for a large-a fixed point in Fig(6j deter- 
mine the fixed points. Closer look shows that RSB fix-points 
lie dense and lead to a fixed point function, well approximated 
by the model function q(a) (dashed grey). 

space can be given by the linear relation among the labels 

777 

{l + k,K= — l + k-1} (4) 

n 

(viewing I > l m i n = Iq = n as the running index) with 
steps of Al = n and m, n, k integer- valued. The choice of 
m/n selects one fixed point of the RSB-flow as k — > oo 
with I — > oo. Steps of Al = n are required to generate in- 
teger values for n (otherwise we wouldn't have numerical 
data). The integer k distinguishes different lines which 
all meet in the same fixed point. Thus the fixed point 
(a*, q*) is labeled by m and n or just by the rational 
number m/n. We have evaluated more than 50 fixed 
points belonging to the exact order function q(a). The 
higher n the larger must be the steps Al, hence one needs 
higher orders of RSB to find enough data points for rea- 
sonable curve- fitting through these points. This is one 
limitation of the method, but the almost linear charac- 
ter of a large number of these lines allows to calculate in 
principle a number of fixed points much larger than the 
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FIG. 5: 



Pade approximants, matching the RSB-flow of dis- 
crete g;, K (a; >re )-data (1 < k < 200) along k — Zl — k/6 — 1 for 
integers k = —4, —3, 4, with initial values lo — 6 + k, kq — 
6 + k and Al — 6. All lines join in the RSB fixed point 
P7,r> = (0.3820,0.2976) (red dot) as « — > oo. (Pade curves are 
displayed here without termination at the fix point. Colors dis- 
tinguish curves, eg red curves: k = 0) 
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FIG. 6: Fixed point example in the a 3> 1-regime is displayed 
with nonlinear RSB-flow along k — 61 + k, k = —2, — 1, .., 2, 
modeled by Pade-curves ending in a large-a fixed point. Its 
position (big blue dot) at {10.6736, 0.9964} and that of a set 
of 50 other fixed points (red dots) , along the model function 
q(a) (green line), is also shown. 



the continuum, while a magnetic field has a similar ef- 
fect on the discrete levels at a = 0. The ratios assume 
the value 1 then. The discrete spectra therefore empha- 
size the critical nature of the points a = and a = oo. 
We present in the following subsections new results for 
these T = levels of parameter ratios and, in Section fVl 
describe their singular finite T crossover. 



1. Level distribution at CP2 

At the critical point CV2 the sub-class of small self- 
consistent parameters eft and a^, which vanish in the 
ooRSB limit (and condense into CV2), obey 



m+2 



21 - 1 
2Z + 1 



and 



H+l 



I 



i + r 



(5) 



with I = k — I and / = 1,2,...; thus the ratios of these 
parameters are discrete and almost equidistant^. Re- 
curring these relations to the smallest parameters of each 
RSB-ordcr k, hence to q K +i and a K respectively, we ob- 
tain 



q K +i-i = (21 + 1)<7k+i, a K -i = (l+ l)a K - 



(6) 



The RSB flow of numerical data up to 200RSB allow to 
conclude that these minimal parameters vanish like 



Qk+1 



1.03059 1.31705 



2.77275 



K 

3.54347 



+ 0(l/n 3 ), (7) 
0(1A 3 )- (8) 



The discretized slope of the order function in the point 
a = 0, assumes the 200RSB value 



1l - Qi-i 2g K+1 



= 0.74345 



(9) 



H-l 



or, by Pade approximation of the RSB flow and extrap- 
olation to oo— RSB, one obtains 



q'(0) = 2 lim K - 



Qk+i 



0.743368. 



(10) 



order of RSB. 



B. Discrete spectra in the k — oo RSB-limit at zero 
temperature 

While the fixed point functions can be derived for all 
pseudo-time values l/o, the points a = and a = oo re- 
main special limits. In a recent article^ we have shown 
that infinitely large subclasses of certain self-consistent 
parameters ratios remain discrete at T = or H = 
even in the continuum limit. These discrete levels reside 
in the limits a = and a = oo when k — oo. Finite 
temperatures lift the discrete spectrum at a = oo into 



As the calculation of fixed points of q(a) in the linear 
small-a regime shows, this agrees with the slope of the 
continuous q(a) for a — > 0. The slope of the order func- 
tion in CV2 provides one almost exact constraint for the 
order function 



q'modei( a = 0) 



1 



vMO) 



1 F 1 a, 7 



e 



w(0) 



0.743368. 



2. Level distribution at CP1 

In the large a limit the characteristic feature are dis- 
crete spectra of 1 — ^-ratios, which are shown in Figure 
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[71 In addition Figure [8] shows that the -i? coefficient of 
the almost continuous order function converges towards 
0.41 except for the largest a-levels. At zero temperature 
the order function differs from 1 by 0.41/a 2 . Thus, ac- 
cording to the large-a expansion of our analytical model, 
the expansion coefficient is constrained to satisfy 



q(a) = l-^!l + 0(l/a 4 ) 
7 a z 

= 1 -0.41^ + O(l/a 4 ), 



(11) 



putting a constraint on a^j. Further constraints can 
be found from very precise numerical characteristics; it 
is planned to use this analysis to narrow down the choice 
of an analytical order function model. 

The discrete spectrum yields a coefficient which differs 
notably from this value, as one can see from Figure [7J 
(right) for the leading divergent ai parameters. 
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100 



FIG. 7: The finite and discrete ratios of 1 — g;-factors, which 
vanish for k — > oo at T = 0, as obtained by Pade extrapolation 
to the RSB limit, plotted versus the level number I. 



C. Approach of equilibrium at T = 0: leading and 
sub-leading scaling contributions 

The nonequilibrium susceptibility Xi is a characteristic 
quantity measuring the distance from the equilibrium so- 
lution at k = oo. The entropy had been seer*2£ to vanish 
like the square of xi- The numerical solutions^ for 
evaluated for all 200 leading RSB-orders, are well fitted 
by the T = 0-form 



xi(/t,r = o) 



0.86 



1.85 



( Ko + K )5/3 ( Ko + K y 



(12) 



0.86«T 5/3 - 1.83k~ 8/3 + 3.12«T 11/3 + 1.85k" 



with kq = 1.278. As discussed in Ref|30j the numerical 
uncertainty of O(10~ 6 ) in the exponent is so small that 



10 20 30 

level number t of the 40 largest a-levels 

FIG. 8: The Figure shows convergence of the large-a behavior 
of a 2 (q(a) — 1) towards ~ —0.41. Some levels of the discrete 
spectrum at a — oo show a deviation. 



the expectation of a rational- valued exponent due to one- 
dimensionality leads to the firm prediction of xi ~ k -5 / 3 . 
The quality and density of the numerical results is even 
high enough to predict the subleading correction and the 
amplitudes as well. 



V. FINITE TEMPERATURE SCALING NEAR 
THE CRITICAL POINT CP1 (a = oo, T = 0, H = 0) 

Naturally one would like to start with a scaling hy- 
pothesis for the free energy F. However the SK-modcl 
has two critical points at T = and the free energy picks 
up contributions from both; in the RSB-limit, it can be 
expressed by integrals over entire crossover range from 
a = (CP 2) to a = 00 (CPi) involving the order func- 
tion q(a). 

Thus it turns out useful to start with the scaling be- 
havior of the self-consistent parameters a; and qi, which 
teaches us how to embed scaling into the order function 
q(a, T) mediating the crossover between the two critical 
regimes. Finally, by expressing free energy and internal 
energy in terms of the order function, and by linking the 
entropy with the non-equilibrium susceptibility, we shall 
arrive at consistent scaling predictions for F, U, and S 
b+elow. 

Let us begin with temperature-normalized block size 
parameters 



ai(n,T) 



mi(K,T) 
T 



(13) 



where we consider first scaling in the (k, T)-plane for 
fixed label /. Wc must analyze the singular behavior 
near the critical point CPI, where diverging a;(«;,T = 
0) — > 00 for k — ► 00 lead to discretely spaced ratios 
uz(oo, 0)/u;_i(oo, 0) in the 00— RSB limit. We identified 
the large order power law divergence 



cli(k, T = 0) ~ k 



5/3 



00 



(14) 
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< a(T) « I/T 

FIG. 9: Main Figure shows the RSB-flow of {ai («), qi(n)} (above) and {ai(ft), (below), for each of the six highest 

parameter levels (I = 1...6) at fixed temperature T = 0.03. Pade approximations (shown for largest level (orange), 2nd-largest 
(green) etcetera) model the flow k —> oo by extrapolation, k — oo is reached in the level crossing point LCV = {au(K) ~ 
18.226, gfc(/t) ~ 0.9986} which separates plateau-regime LCV < a < 1/T from the oo-RSB continuum < a < LCV. The Parisi 
box size m(T = 0.03) = 0.5467 is the break point value. Inserts illustrate that qi- and a;-levels meet in the LCV for k = oo. 



for the subclass of large parameters a; (their number also 
grows to infinity as k — > oo). 

The linear temperature decay of all Parisi box sizes 
mi(K,T) ~ T holds for all finite RSB-orders, but not all 
m's should vanish in the RSB limit at zero temperature, 
since the break point is not expected to be at mi = 
(even in the T -> 0-limit^). Thus, one should describe 
a non-commuting limits T — ► and K — > oo properly. 

The Taylor series, valid as a low temperature expansion 
for any fixed finite RSB-order, 

mj(«, T) = oj(«, T) T = 0)T + ^(s, 0)T 2 + 0(T 3 ) 

(15) 

will anyway break down for those levels Z for which the 
expansion coefficients diverge as k — > oo. In accordance 
with the anomalous power law (|14[) it will be shown below 
by means of the fixed point order function that the cor- 
rect scaling form for this CPi-divergent parameter sub- 
class reads 



50-RSB and satisfies 



a l ( K ,T) = ^ 3 f ai (T/ K - 5 / 3 ) 1 



(16) 



where the scaling function is well approximated by a low 
order (2, 3)) Pade series (one may also use (1, 2) or (3, 4) 
series) 



hi 0) = 



coj + CllX 



(17) 



1 + dijx + d,2,ix 2 
This form fits well the available finite T data up to 



1 



/a,(0) = Co.; finite and f a ,(x) ~ - for x — > oo. 

x 

(18) 

The crossover line can be described by the character- 
istic (crossover) temperature 



(19) 



Beyond the crossover line, for temperatures T 3> 
Ti(k), the box sizes mi(x) = x f ai (x), which belong 
to the CPl-divergent sub-class of et/'s, approach finite 
temperature-independent values. One obtains 



lim mi(x) = cij/ckj. 



(20) 



While direct fits of our numerical data yield already a 
crude estimation of mi(oo) for the break point, it was 
mentioned in Refl30l that 50-RSB is not sufficient to de- 
termine the break point for arbitrary low temperatures. 
Yet, for T = 0.015 a reliable break point value was de- 
termined by another procedure. 

Here we arc interested to obtain a good estimation 
of the breakpoint in close connection with the scaling 
picture. Therefore we employ the fixed point method and 
indeed succeed in finding a good approximation down to 
even lower temperatures and also answer the question 
whether the limit m/(oo) in Eq. (|20p shows a level index 
dependence or not. For finite mi an ^-dependence would 
have implied a discrete distribution. We shall find in 
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subscctio riV Al that all ratios become level-independent 
in the large x limit 

rrnW adx) / a ,M for 
mi-i(x) ai-i{x) Jai-Ax) 

(21) 

The crossover from discrete parameter spectra for T <C 
7i(re) to the continuum on the other side of the crossover 
line, for T 3> Ti(re), is a rather singular effect mediated 
by the scaling function. We introduced above a scaling 
function which allows to suppress the discrete spacing 
between q- and a-parameters as one moves through the 
crossover line 2i(re) ~ re~ 5 / 3 . 



A. Forbidden level crossing at finite temperatures 
determines the break point 

We employ now the RSB-fixed-point technique to ex- 
tract approximate values for the break point for rather 
low temperatures. 

For this purpose, we consider fixed finite temperatures 
T and fixed level numbers I (down to lowest T and I small 
to catch the diverging-a subclass near CV1) and study the 
RSB-flow of the solutions {aj(re), gz(re)}, and also those 
of the complementary type {a; (re), g/+i(re)}, from low or- 
ders up to re = 50 as illustrated by Fig|5]for an arbitrarily 
picked temperature T = 0.03. Pade-approximants fit the 
RSB-flow well and these extrapolated curves meet pre- 
cisely in the same point. These cu+rves would cross each^ 
other, but then violate the reality condition of the self- 'yf 
consistent method beyond the level crossing point. We *J 
consider the level crossing point therefore as the limit of ^ 
the nontrivial part of the order function, hence as the 
breakpoint. 

The scenario remains the same for arbitrary fixed tem- 
peratures, only the extrapolation range increases with 
the level number and therefore becomes less accurate 
for smaller temperatures. Yet reliable solutions were ob- 
tained down to temperatures T w 0.005. The Figure In- 
sert emphasizes the fact that the solutions indeed reach 
the level crossing point as re — > oo. 

Approaching zero temperature and the RSB limit 
along the crossover line, x fixed, with Ti(re) ~ re~ 5 / 3 , 
leads to a discrete set of different Parisi box sizes mi(n — 
oo,T = 0). 



B. Nonequilibrium susceptibility \i 

The scaling form of the non-equilibrium susceptibility 
Xi( K i T) can be given in terms of a scaling function f\ by 

X i(n,T) = K- 5 / 3 f 1 (T/ K - 5 / 3 ) (22) 

where fi(x) ~ x, x — > oo, and /i(0) 53 0.86, reproduces 
the data and the leading re-decay at T = 0, as in Eq. (fT2")) . 



VI. MAGNETIC FIELD SCALING AT 
CRITICAL POINT CV2 (DIVERGING 
PSEUDOTIMES l/o -» oo) 

The magnetic field dependence at T = is expected to 
yield a plateau-like cutoff of the order function of simi- 
lar shape as described in the Parisi form q{x). We study 
now the field dependence of the smallest order parameter 
q K+1 (H,T = 0) in re-th order of RSB. 20 orders of RSB 
turn out to be enough to extract the exponent describing 
the decay of q K +i as the order of RSB tends to infin- 
ity. Guided by the results of finite temperature, where 
one single non-trivial rational exponent appeared, we ob- 
serve an exponent 2/3 to provide a reasonable picture for 
extrapolation towards oo RSB. 

We first identify the qi ~ re _1 -law for (infinitely many) 
order parameters which vanish as re — > oo. The scaling 
hypothesis for (re, H )-scaling, valid for the vanishing or- 
der parameters (ft (re, H), can be formulated as 

1 /H 2 / 3 \ 
q i ( K ,H,T = 0) = -f i (-—) (23) 
re \ 1/re / 

with /i(0) ^ and fi(x — » oo) ~ x. 
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FIG. 10: Upper and lower bounds of the order functions 
are shown for {20RSB,H = 0.3} and in 200RSB for zero 
field. For comparison the H = fixed point order function 
(green) is included. The discrete 20RSB (black dots) order 
(interpolation-)function and two fixed points obtained from 
the 1 — 20 RSB-flow extrapolated to re = oo (light-grey lines) 
in the finite field H = 0.3 are displayed in addition. The fixed 
point at a*pi — 0.568, g^; = 0.448 defines the plateau-cutoff of 
q(a) in the finite field H = 0.3. 

The numerical procedure chosen in order to arrive 
at this proposal has been to extrapolate to ooRSB the 
smallest gi-values at fixed non- vanishing small magnetic 
fields. The higher the field the less orders of RSB are 
needed (similar as in the case of finite temperatures). 
Twenty steps of RSB generate almost exact results down 
to H w 0.15. Extrapolation of the RSB-flow is hence re- 
liable down to much smaller field- values, where one has 
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FIG. 11: Numerical data for the magnetic field range < H < 
0.5 are shown from 1st to 20th order of RSB (green). Dots 
(dark red) show the plateau-height small-q cutoff obtained 
from the fixed point order function. At H = dots show the 
calculated orders up to 200RSB for comparison. The RSB 
flow for q m in{K, H — 0) is given by a 1/k law, while q m in(n = 
oo,H) obeys an i/ 2//3 -law. 



already entered the critical regime. Thus many RSB fixed 
point values (at k = oo) are well approximated and can 
be used to match a power law w.r.t. the magnetic field. 
In this way the magnetic field exponent of Eq. ([23]) is 
found to differ only by 0.003 from the value 2/3 which 
led to the assumption that this rational number is exact. 



VII. SCALING BEHAVIOR OF THE FREE 
ENERGY F, INTERNAL ENERGY U, AND 
ENTROPY S 

Low temperature expansions of internal energy U, en- 
tropy S, and the free energy F = U + TS were reported 
in the framework of our high order RSB analysis, and 
found in agreement with already known results. In the 
present context of scaling theory, we also look for scaling 
of RSB-parameters together with temperature and also 
small field variation. A useful way to study the RSB-flow 
in terms of k— scaling is by invoking the internal energy 
formula at T = and H = 

U(k, T = H = Q) = - X x-\Y j ai (qf - qf +1 ) (24) 



i=i 



lim U(k,0) 



1 

2 Jo 



da (1 - g(a) 2 ) 



(25) 



The summation includes contributions from both crit- 
ical points and from the crossover regime in between. 
Consequently one cannot expect to obtain scaling laws 
from a single hypothesis imposed on the total free en- 
ergy. The problem has more in common with critical 
dynamics, however with two critical points in the long 



pseudo-time limit (1/a — » oo) and in the short pseudo- 
time limit (1/a — > 0). 

As reported in Refl30l the free energy has a low T ex- 
pansion in the RSB limit given by F = F(T = 0)-S(T = 
0)T + J2k=2 fkT k , where the leading temperature behav- 
ior is F(k = oo, T) - F(oo, 0) ~ T 3 . The leading large-K 
correction of the T = free energy has also been reported 
to decay like k -4 . 

In the large-a regime, temperatures scale like k -5 / 3 
and hence the large-a scaling contribution is SF ~ re~ 5 . 
Thus the leading temperature dependence belongs to a 
sub-leading SF ~ k~ 5 correction. 

We attempt to distinguish singular scaling contribu- 
tions from both critical points from the non-singular con- 
tributions to the free energy. The small a-regime contri- 
bution can be estimated from 

F(k = oo, T = 0) = U(k = oo, T = 0) = E {H) 

i r°° 

= -- da (1 - q 2 (a)) - M(H)H, (26) 

where M(H) denotes the field-generated magnetization. 
Recalling the small-a expansion of the order function, 
q(a) ~ a — const. a 3 + 0(a 5 ), one must expect an iJ 10 / 3 - 
contribution from the plateau-regime, which implies also 
an 0(k~ 5 ) contribution. The free energy data are com- 
patible with an iJ 10 / 3 small field scaling part. 

It must be concluded that the leading correction k~ 4 
must originate in the intermediate a- regime (not yet 
identified in detail). It can, after all what was said before, 
not be assumed to be a scaling contribution. We should 
therefore attribute it to the regular free energy part. 

The entropy was found to obey^i 

S(k, T = H = 0) = -- A xi (k, T = H = 0) 2 . 

It is known that only the large-a regime near CV1 is re- 
sponsible for the leading K-behavior of xi at zero temper- 
ature, hence this holds also for the T = 0-cntropy. Since 
thermal behavior is also caused by the CV1 contributions, 
we can therefore claim that the scaling-contribution to 
the entropy obeys 

^( K ,T) = K - 10 / 3 / s (T 2 / K - 10 / 3 ), (27) 

with fs(x) = —0.72 x (see Refl30h for x — > oo and 
fs(x) ~ —0.185 for x — > 0. Thus the entropy contributes 
to the leading 0(T 3 ) low temperature correction of the 
free energy. This TS'-term contributes again only a sub- 
leading correction SF ~ k~ 5 from the large K-scaling 
regime. 

Let us recall the large-re dependence of the free energy 
at zero temperature, well described by the optimal fitting 
form 



F(k,T= 0) = F(oo,0) 



c 4 



f'5 



(k + k ) 4 {k + k ) 5 



where excellent Pade-fits yield the constant kq = 1.28. 
The leading correction k does neither originate from 
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the scaling regime near CV1 nor from that near CV2, and 
hence must be expected not to scale. We therefore con- 
sider it as part of a regular ^-contribution F reg (n, T, H). 

Thus we propose that the free energy consists of a sum 
of a regular and of two singular parts, where the latter 
ones scale according to whether they are CV1- or CV2- 
critical. 

As a consequence of this two-critical point picture and 
in agreement with the numerical data, we separate two 
singular contributions, which offer different scaling be- 
havior, from a regular part F reg by 

F(K,T,H) = F reg ( K ,H,T) + Ff^) ( K ,T) + F^ (k,H) 

(28) 

where the magnetic-field controlled critical point CV2 
and the temperature-controlled critical point CV1 con- 
tribute respectively 

F^(k,T) = K - 5 f cpl (T/ K - 5 / 3 ), (29) 

and 

F^(k,H) = «- 5 / cp2 (^ 2/3 /«- 1 ) (30) 

with f cp i(x) ~ x 3 , f cp 2 ~ x 5 for x — > oo and both fi- 
nite for x — > 0. This claim refers to the leading scaling 
behavior at CV1 and CV2; corrections with analytic In- 
dependence near CV2 and analytic field-dependent cor- 
rections near CV1 may occur. 

A contribution — ^x( K )H 2 -term, which yields the lin- 
ear equilibrium susceptibility from —&fjF, belongs to the 
regular part F reg with x( K — > oo, T < T c ) = 1. 

It is interesting trying to translate the given power laws 
into scaling with the number N of spins for the finite N 
SK-models^ which corresponds to a finite size system 
with N = L d , d denoting the real space dimension. Scal- 
ing with L or N delivered a leading correction ~ iV~ 2 / 3 
for the finite SK-mode l 17 ' 18 . If we would assume scaling 
of the leading correction k with N, a scaling function 
depending on N^ 1 ^ 6 / kT 1 would resul t 37 ' 38 . However this 
rests on the assumption that the leading N~ 2 / 3 energy 
correction arises from the entire a-rcgime. Many open 
questions seem to show up here. 

VIII. FIXED POINT DISTRIBUTIONS 

A. Ground state energy Eq 

We can extract more detailed information from our nu- 
merical analysis of RSB in the SK-mode l 20 ' 30 beyond the 
calculation of the global ground state energy. The RSB- 
flow of the energy level distribution and naturally the 
energy density eo(a) as a function of pseudo-times can 
be given. In the latter case, a test of our analytic order 
function model against the numerical results 3 ^ is pro- 
vided by the use of q(a) and of q'(a). Both are required 
in the ground state energy formula in Eql26l according to 

/>oo />oc 

Eq = da eo(a) = — / da aq'(a)q(a). (31) 
Jo Jo 



Using the analytic form (3) and high RSB-order results 
for k = 100, 110, 120, ...200, we obtain FiglH^ 7 . We do 
not find exponential tails in this energy distribution, in- 
stead observe simple power law decay in the limits of 
small and large a. 

A second important representation shows the energy 
level contributions from 

e (l, n) = -- Urn ai(K) (qf (k) - , (32) 

as a function of normalized level index 1/k, and with 
boundary conditions do = P, qo = 1- The sum over 
all energy levels eo(l,n) with level index / = 0, 1,2, ...k 
for each calculated RSB-order yields the RSB-fiow of the 
ground state energy 

K 

£*,(«) =X>(i,K) (33) 

1=0 

towards the exact valued E (k = oo) = Eq. A proper 
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FIG. 12: The (negative) energy density — eo(a) of the ground 
state energy Eo = J °° da eo(o), obtained from our analytical 
model function q(a) as — e(o) = a q'(a)q(a) (black underlying 
curve), is shown to agree well with discrete numerical 200- 
RSB results . 

normalization of level numbers by the RSB-order k, dis- 
plays the level-distributions for each RSB-order on the 
same interval of unit length. Subsequent rescaling of the 
energy level allows to visualize the RSB-fiow towards one 
fixed point energy distribution (which of course depends 
on the rescaling factor—) . FigfT3] shows two choices (l- 
and k- rescaling of eo(l,K)) - in both cases the convergence 
towards the fixed point function is obvious. 

Fixed points (under RSB-fiow) have been calculated 
in the same way as shown before for the order function. 
For example, fixing 1/k to a rational number m/n within 
the unit interval, one can see many of the leading fixed 
points in Fig[]2] following the RSB-fiow along vertical 
lines fixed by m/n. The piecewise dense set of calculated 
fixed points was obtained by an extrapolated Pade ap- 
proximation for n — 2, 51 with m = 1, n— 1. These 
fixed points are shown in Fig ll3l together with their fit 
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FIG. 13: Part (a) shows ground state energy solutions (dots) 
for k eo(l, k) (or I 6o(1,k)) and their RSB-flow from k — 10 to 
k = 200, Ak = 10, towards their fixed point functions [C]p* (C) 
on the unit interval of levels 1/k — > C- 00 shows differences 
Ai,2 between 200RSB (interpolation functions [ClPehoo) and 
the fixed point functions [C]Pe(0 (red curves in (a)). 

function, obtained here as an (8, 8)-Pade series. The fixed 
points are piecewise dense with some gaps near 'lead- 
ing' fixed points (eg at , which become however closed 
as one higher orders are evaluated. The fit function (in- 
terpolating between the dense regions) represents an ap- 
proximation for the exact fixed point energy distribution 
function £q(£) with 1/k — > £ in the oo-RSB limit. The 
numerical integration of the approximated function p* (£) 
(which corresponds to e*(a) of Eq. (|3"Tj) transformed form 
< a < oo onto the unit interval < C < 1) yields^ 

E* = E*(T = 0) = / d( p* e (0\approx « -0.76314. 
Jo 

(34) 

By reproducing the correct valued up to O(10 5 ), this 
provides a good test of the fixed point method. An 
alternative calculation, using Eg. pip with plugged in 
fixed point order function confirms the numerical value 
Eq. The inserted figure shows the magnitude of energy- 
corrections per level I occurring from 200-RSB to the 
exact oo-RSB energy per level (recall that / labels the 
Parisi boxes of the RSB order parameter). 

Different power law decays arc observed in the small 
1/k (CV1) and in the 1/k w 1 range near (CP2). 

An analytical modeling of the fixed point energy dis- 
tribution must be attempted in the future; it might 



reveal more valuable information about the relation 
with directed polymers and/or with the KPZ-universality 
classes^ 4 -. 

Energy distribution functions play an important role 
in the characterization of directed polymer a 40 ' 41 i 42 . Gen- 
eralized Gumbel statistics 4 ^ were found to describe the 
statistical fluctuations of global quantities (like the en- 
ergy). It is perhaps in this respect where a clear dis- 
tinction between the directed polymers and the present 
universality class can be made. But this detailed compar- 
ison is beyond the scope of the present paper and should 
be attempted in the future. 

B. Equilibrium susceptibility per level 

To conclude this section we extend the described 
method to the x( a )~density of the equilibrium suscep- 
tibility x an d in particular to the distribution per level 
I. In the RSB- limit, the total x is known to be equal to 
1 in the entire ordered phase. The RSB flow thus moves 
towards a fixed point function x( a ) = a <l'( a ) with the 
property J Q dax(a) = 1 (this had been used before as a 
constraint for our analytical order function mode l 20 i 27 ). 

Let us now study the RSB-flow of the discrete repre- 
sentation x(Z,k) = ai (k) (qi (k) — qi+i(K)). The result 
for the susceptibility per level / (normalized by RSB- 
order k) is shown in Fig 1 141 (and corresponds to the en- 
ergy per level distribution shown in the preceding fig- 
ure). The shape recalls universal distributions observed 
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FIG. 14: RSB-flow of the equilibrium susceptibility per level 
I from k — 10- to k — 200 (in steps of 10 orders of RSB and 
rescaled by I — > 1/k to a unit interval) towards the fixed point 
susceptibility distribution p£(C) (black curve) 

for the KPZ growth processes 4 ^. This relation or map- 
ping must be studied in the future, particularly because - 
as explained in RefH - the related statistical fluctuations 
have been associated with universal critical behavior. 

Beyond the flow of the finite RSB orders K = 
10, 20, 30, 200 we have added the fixed point function 
pt{() for the susceptibility density (ic x( a ) transformed 
onto the unit interval < £ < 1) which must obey 
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J dC, = 1- A simple approximate calculation of 

the interpolating fixed point function reproduces the ex- 
act constraint with an error of only O(10 -5 ). Again this 
confirms the power of the method, which can eg be used 
to test analytical proposals. 

Small changes from 200-RSB to oo-RSB are resolved 
in FigfTJ] and in FigfT^l Their tendency is to make the 
distribution more symmetric. Yet the distribution per 
normalized level remains asymmetric as for the energy 
distribution (as a function of the dense levels l/n), a uni- 
versal fact that has been observed as a special feature 
of the SK-model in contrast to symmetrical distributions 
finite-range spin glasses. 



IX. SCALING WITH THE 
PSEUDO-DYNAMICAL VARIABLE OF THE 
ORDER FUNCTION q(a) 

In previous publications we found a Langevin-type 
representatio n 20 ! 33 for a logarithmic derivative of the or- 
der function q(a) with respect to 1/a. This ordinary dif- 
ferential equation (without stochastic field) is much sim- 
pler than the exact partial differential equations, which 
is a consequence of the existence of scaling behavior and 
of homogeneous functions. It is well-known that scale 
invariance and the so-called similarity method reduces 
partial to ordinary differential equations^. Therefore, 
at least near the critical points one can expect ordinary 
differential equations to describe RSB. 

The Langevin-type of differential equation could how- 
ever be reshaped in terms of different pseudo-dynamic 
variable a, l/o or other forms. The differential equa- 
tion remains to be relaxational and thus there remains 
some arbitrariness in the choice of the proper 'time' vari- 
able r. If we wish to apply dynamic scaling to the RSB- 
representation, we are unfortunately bound to make a 
definite choice. Let us consider o + l/aasa pseudo-time 
in order to conform with the expectation that critical 
behavior at either of the points CVl or CV2 should oc- 
cur in the long-time limit. Then at CVl we would get 
t ~ a — > 00 while r ~ - — > 00 at CV2. 

a 

We may now consider pseudo-dynamical scaling by 
studying the a-dependent quantities like the order func- 
tion near CVl and CV2. 

Near CVl the order function obeys 

q(a, k,T) = 1 + a- 2 f q (T 2 /a- 2 , a 2 /« 10/3 ) (35) 

with f q (x,0) ~ x, f q (0,x) ~ x, and f q (0, 0) fi- 
nite. In terms of the transformed order function (ft ~ 
d 1 / a log{q{a)) one gets (ft ~ 1/a at re = 00, T = and 
~ T at a = 00, k = cxd. This would allow to extract an 
exponent /3 = 1, and together with £ K ~ a 3 / 5 ~ T~ 3 / 5 
the correlation exponent v = 3/5 results. Then, using 
r~ 1/aasa pseudo-time variable near CVl, the dynamic 
exponent follows from r ~ as z = 5/3, remarking also 
that z v = 1. 



Given the already mentioned similarities with directed 
polymers, the known relationship between those and the 
KPZ-equation^i suggests a comparison between pseudo- 
dynamics of RSB in the SK-model and the dynamic KPZ- 
behavior. 

We note that dynamic critical exponents were recently 
reported by Canet and Moored for two universality 
classes of the KPZ-equation. For one type of approx- 
imate solution of the Flory-Imry-Ma or RSB-type, the 
dynamic exponent assumed the value z = (4 + d)/3 be- 
low two dimensions, hence z = 5/3 in d = 1. Hence 
we state that the exact pseudo-dynamic critical expo- 
nent of RSB in the SK-model maps to the one of FIM 
or RSB-approximate solution of the KPZ-equation in ID 
(provided one agrees to make the choice of 1/a being 
the pseudo-time which corresponds to the real time of 
KPZ). As in the DP-analogy, this should correspond to 
the domain- wall solution and hence to CVl. 

There is however also the known exact result of the ID 
KPZ-equation z = 3/2 also given by Canet and Moored. 

One may suspect that this result should be mappable 
to pseudo-dynamic behavior of RSB-SK near the second 
critical point CV2. Indeed, if we would conserve z v = 1, 
the same exponent z = 3/2 would be obtained near CV2. 
We do not have any reason for this choice, and the explicit 
scaling of the order function near CV2 does not confirm 
this value, neither for the choice r = 1/a nor for r = a. 
This question must remain open. 



X. DETAILED STRUCTURE OF THE ORDER 
FUNCTION DERIVATIVES q'(a) AND q"(a). 

The derivatives depend much more specifically on the 
pseudo-time variable than q(a) itself. Failure of an an- 
alytic model function becomes detectable more easily in 
the derivatives. In order to control our modeling, we 
studied analytical fits first of all 200-RSB data, and sec- 
ondly of the 50 calculated fixed points. Taking q'(a) di- 
rectly form the analytical form q(a) as given by Eq.Q, 
we find good agreement with the discretized slope cal- 
culated from the fixed points. This is demonstrated in 
the main part of Fig|T5j In addition the insert shows the 
second derivative d 2 q(a), where the two analytic mod- 
els (red and blue curves) show a small difference. We 
note in passing that the shape of the 2nd derivative q" (a) 
shows a similarity with the 2-loop correction in the Y- 
correlator of (1 + l)-dimensional random bond pinned 
manifolds^ (we don't know whether this similarity has 
a deeper reason). The maximum seen in q'(a) expresses 
the Crisanti-Rizzo curvatur o 20 ' 50 , a slight non-linearity 
of the order function in the small a regime. It is however 
this contribution, which renders an analytical fit rather 
awkward. An analytical model which fits well the neigh- 
borhood of the critical points a = and a = 00 can 
have a simpler shaped, but we want to get the pseudo- 
dynamic crossover right as well. Global quantities like 
the energy (integral over all a), picking up only small 
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FIG. 15: The characteristic maximum in q'(a) is analytically 
well reproduced. Red dots show the discretized derivative 
obtained from 50 fixed point {a*,q*}. The insert shows the 
second derivative of the model order function q(a). 

contributions nearby the critical points, depend on the 
crossover regime modeling. This can be seen in Eq. (j3"Tj) 
as well as in Fig. 12 for the energy density. 

XI. CONCLUSIONS 

In this article we formulated a scaling theory of the 
flow towards full replica symmetry breaking (RSB) at 
T = 0, for finite temperatures, and for finite magnetic 
fields in the SK-model. Several fixed point functions of 
RSB-fiow were evaluated. 

The analysis was guided by 

1. a large set of high-precision numerical data, with up 
to 200 sclf-consistently solved orders of replica symmetry 
breaking for the T = SK-spin glass and still a high 
number of orders for finite temperatures and magnetic 
fields, 

2. by the identification of two critical points (at zero 
temperature and zero magnetic field), which are distin- 
guished by two different pseudo-dynamic limits as ob- 



tained in an analytic picture of a Langevin-type equation 
in ReflHl! and 

3. by representing nonanalytic behavior near each of 
these critical points in the framework of the scaling the- 
ory of critical phenomena. 

Power laws and scaling functions were identified by 
fitting the leading 200 RSB-orders of self-consistent so- 
lutions deep inside the SK spin glass phase; non-integer 
exponents were found and identified as rational numbers, 
characteristic of one-dimensional RSB-bchavior. This 
ID-character originates in correlations on the pseudo- 
lattice of RSB-orders k. By means of scaling functions 
we demonstrated how these nonanalytic ID-correlations 
enter in temperature- and field-dependent power laws in 
the ordered phase. 

The universality class of replica symmetry breaking in 
the SK-model called for comparison with other physical 
systems, and shows similarities with directed polymers. 

The decoupling of a magnetic field sensitive critical 
point from a temperature-sensitive one was embedded in 
an unconventional scaling hypothesis for the free energy 
and found to be consistent with the numerical data. 

The RSB flow was used to generate an order parame- 
ter fixed point function, serving as a crossover between 
the two different pseudo-dynamical critical limits. Its 
fine structure was revealed by the leading derivatives, 
again confirming excellent agreement between analytical 
model and fixed point function. 
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